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ABSTRACT 


Turbularization  of  an  acoustic  boundary-layer  (Stokes  layer)  on 
impermeable  and  permeable  surfaces  is  analytically  considered.  The  theoretical 
approach  utilizes  a  second-order  closure  model  of  turbulence.  Both  an 
approximate,  closed-form  solution  and  a  more  comprehensive  finite  difference 
solution  of  the  time  dependent,  parabolic,  one-dimensional  governing  equations 
are  obtedned.  For  simple  acoustic  boundary-layers  on  impermeable  surfaces,  both 
the  approximate  solution  and  the  numerical  results  for  the  critical  acoustic  Mach 
number  required  for  turbulent  transition  are  qualitatively  confirmed  by 
experiment.  The  calculations  for  acoustic  boundary-layers  with  transpiration 
(injection)  indicate  a  substantial  reduction  of  the  acoustic  Mach  number  required 
for  transition,  up  to  a  limiting  injection  velocity  that  is  frequency  dependent.  The 
results  may  provide  a  mechanism  for  flow-related  combustion  instability  in 
practical  systems,  particularly  solid  propellant  rockets,  since  turbularization  of 
the  near-surface  combustion  zone  coiild  result  at  relatively  low  acoustic  Mach 
numbers. 

An  analysis  of  the  transitional  and  turbulent  reactive  acoustic  boundary 
layer  on  a  homogeneous  solid-propellant  surface  is  conducted  to  investigate 
potential  mechanisms  of  combustion  instability.  A  new  technique  is  developed  for 
the  condensed-phase  thermal  layer,  in  which  the  propellant  space  is  mapped  onto 
the  gas  space  and  efficiently  solved  using  the  same  adaptive  numerical  grid. 
Results  are  obtained  at  an  acoustic  pressure  node  in  the  absence  of  a  mean  axial 
flow.  The  results  indicate  that  acoustically  induced  transition  can  occur  at 
relatively  low  acoustic  pressure  amplitudes,  propellant  response  to  harmonic 
axial  velocity  fluctuations  is  both  rectified  and  displays  a  mean  augmentation 
("D-C  shift"),  and  that  nominal  propellant  combustion  parameters  lead  to 
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increasing  susceptibility  to  acoustic  transition  at  elevated  mean  pressures. 

Further  analysis  of  the  the  turbulent  reactive  acoustic  boundary  layer  on  a 
homogeneous  solid  propellant  surface  is  conducted  to  investigate  the  variation  of 
acoustic  erosion  along  the  length  of  the  chamber  of  a  solid  rocket  motor  or  center 
vented  T-bumer.  An  order  of  magnitude  analysis  indicates  that  the  enthalpy  and 
species  equations  remain  similar  in  the  presence  of  acoustic  pressure 
oscillations.  Results  are  obtained  throughout  the  length  of  the  chamber  in  the 
absence  of  a  mean  axial  flow.  The  results  indicate  that  the  increase  in  mean  burn 
rate  (DC-shift)  observed  in  many  T-burner  experiments  may  be  due  to 
turbularization  of  the  acoustic  boundary-layer  above  the  burning  propellant.  A 
comparison  is  made  between  predicted  and  experimental  response  functions  for 
acoustic  and  steady-state  erosive  burning.  The  threshold  velocity  for  acoustic 
erosion  is  foimd  to  be  significantly  lower  than  that  of  steady-state  erosion,  and  the 
amplitude  and  harmonic  content  of  the  propellant  response  varies  considerably 
throughout  the  chamber. 
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NOMENCLATURE 


ao  imperturbed  sonic  speed 

As  pre-exponential  constant  in  surface  reaction 

Bg  pre-exponential  constant  in  gas  phase  reaction 
Cp  specific  heat  at  constant  pressure 
Ctc  specific  heat  of  condensed  phase 
f  frequency,  Hz 

h  specific  sensible  enthalpy 

h°  heat  of  formation  of  species  a 
H  total  enthalpy,  h  +  u’ui/2 
k  thermal  conductivity 

kg  equivalent  sand  roughness  height 

Lg  heat  of  (dissociative)  evaporation  or  sublimation 
M  Mach  number  defined  using  ao 
n  normal  burning  rate  pressure  exponent 
p  static  pressure 

q  turbulence  intensity,  (u'*Ui)l/2 

r  radial  distance  from  centerline 

r  condensed  phase  surface  regression  rate 

R  inner  radius  of  a  cylindrical  duct 

Reac  axial  acoustic  Reynolds  Number,  u"c/(f/vc)^^2 
Reg  injection  Reynolds  number,  PgVgS/pg 
Ret  turbulence  Reynolds  number,  pqA/p 
Ru  universal  gas  constant 

t  time 

T  static  temperature 

Ta  activation  temperature 
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uj  velocity  vector  (u,v,w) 

Wav  average  moleciilar  weight 
X  j  coordinate  vector 

y  distance  from  surface,  6  -  r,  (planar  flow);  or  R  -  r,  axisymmetric  flow 
Y(x  species  mass  fraction 
a  index  for  chemical  species 

thermal  diffusivity  of  the  condensed  phase,  /( c^^ ) 

Pg  temperature  exponent  of  gas  phase  reaction  pre-exponential  coefficient 
6  a  non-injected,  laminar  acoustic  boundary-layer  thickness,  (v/f)l/2 

5aa  actual  injected  acoustic  boundary-layer  thickness 

6f  normalized  flame  height,  kg/ m  Cp 

Ahg  heat  of  reaction  per  xinit  mass 

condensed  phase  coordinate 
<()  characteristic  length  scale 

O  concentration  (pressure)  exponent  in  reaction  rate 

A  turbulence  macro-length  scale 

n  relative  static  pressure  amplitude  I  p"  I  /  p 

p  viscosity 

V  kinematic  viscosity 

p  density 

a  =k/cp 

Ov  =  [v'8V's/v82]1/2 

X  normalized  time,  f  x  t ,  cycles 

u  index  for  planar  or  axisymmetric  geometry 

0)  specific  reaction  source  term 

^  linear  stretching  factor  used  in  condensed  phase  mapping 
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average  of  variable  over  turbulent  fluctuations 
<>  time  mean  of  variable 

'  turbvilent  fluctuation  value  of  variable 
"  acoustic  fluctuation  value  of  variable 

*  reference  condition 

Subscripts 
a  acoustic 

c  duct  centerline 

e  value  at  edge  of  boundary  layer 
f  value  at  edge  of  flame  zone 

g  gas  phase 

h  condition  at  port  head  end 

1  value  of  variable  at  local  axial  position 
m  maximum  absolute  value 

s  condition  at  surface 

ss  steady-state  value  of  variable 

71  condensed  phase,  propellant 
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SECTION  1:  INTRODUCTION 


1.0  OVERVIEW  -  ACOUSTIC  INSTABILITY 

Oscillatory  flows  in  ducts  can  be  sustained  by  a  variety  of  interactions. 
These  range  from  purely  fluid-dynamically  induced  motions  (e.g.,  large  scale 
vortex  shedding  or  shockwave  instability),  to  the  more  energetic  motions  possible 
with  diabatic  flows.  Particularly  severe  oscillations  can  occur  in  solid  propellant 
rocket  chambers,  for  example,  wherein  the  oscillatory  motion  can  be  driven  by 
interactions  with  the  substantial  energy  release  inherent  in  near-surface 
combustion  processes. 

Prior  analjrtical  work  in  "combustion"  instability  in  solid  rockets  has 
identified  some  of  the  mechanisms  which  can  produce  velocity-field  coupling  to 
the  overall  instability  process  (see,  for  example,  the  review  by  Culicki).  These 
studies  have  shown  that  for  a  simple  longitudinal  standing  acoustic  wave  in  a 
duct,  a  Rayleigh  diabatic  instability  criterion  can  result  which  is  dependent  upon 
the  oscillatory  motion  of  the  gas  column.  Since  the  energy  release  occurs  in  the 
near  surface  region,  analytical  work  has  addressed  acoustic  boundary-layer 
effects  on  propellsmt  combustion  response  (see,  for  example,  Lengele^).  Implicit 
in  several  response  function  analyses  is  the  assumption  that  the  acoustic 
boundary-layers  behave  quasi-steadily  and  in  phase  with  the  longitudinal  acoustic 
velocity  outside  the  boundary  layer. 

Recent  work  on  combustion-flowfield  interactions  in  solid  rocket  chambers 
has  analytically  and  experimentally  examined  fundamental  fluid-dynamic 
aspects  of  mean  flow,  acoustic  wave  and  turbulence  behavior.  Hydrodynamically 
modeling  the  flow  with  a  semi-enclosed,  porous-walled  duct  (with  large  injection 
through  the  wall),  Brown  et  al3.4  experimentally  confirmed  the  transitional 
behavior  of  the  mean  flow  predicted  in  Ref.  5.  The  major  emphasis  of  this  non- 


reactive  flow  experiment  was  to  investigate  the  effect  of  low  and  moderate 
amplitude  forced  acoustic  oscillations  on  the  flow,  using  both  hot-wire 
anemometry  and  surface-mounted  hot-film  sensors. 

Basing  their  conclusions  on  the  surface  sensor  measurements,  Brown,  et 
al.  suggested  that  near  surface  turbulence  produced  by  the  mean- flow  transition 
process  appeared  to  "destroy"  the  coherent  response  of  the  sensors  downstream  of 
transition.  Consequently,  for  low  acoustic  amphtudes,  a  surface  response  capable 
of  inducing  instability  would  most  probably  originate  in  the  head  end  (pre¬ 
transition)  section  of  a  rocket  chamber.  However,  it  is  noted  that  for  large  axial 
distances,  low  amplitude  acoustic  signals  appear  to  be  resurrected  in  the  post 
transition  region.  Lower  amplitude  harmonics  of  the  driving  frequency  were 
observed  by  Brown  in  both  the  surface  and  hot  wire  measurements.  Their 
existence  is  consistent  with  the  nonlinear  behavior  of  a  (potentially  thick)  acoustic 
boundary  layer  (Stokes  layer). 

Analyses  of  laminar  acoustic  boundary-layer  phenomena  involved  in 
propellant  response  have  been  offered  by  Flandro®;  Click  and  Renie’^;  Ben  Reuven® 
and  Hedge,  et  al.9.  Zinn  and  colleaguesiO  have  obtained  Schlieren  and  radiation 
measurements  for  a  premixed  reactive  Stokes  Layer  on  a  porous  plate.  Baum  and 
Levine^ ^  utilized  a  full  unsteady  Navier-Stokes  solution  method;  their  results 
showed  a  Stokes  layer  persisting  even  in  the  presence  of  strong  injection 
velocities.  The  thickness  of  the  Stokes  layer  was  of  the  same  order  as  the  Stokes 
layer  for  noninjected  flow.  The  presence  of  vortices  above  the  acoustic  boundary 
layer  were  also  indicated  in  the  results.  Although  several  aspects  of  these 
analyses  are  interesting,  the  following  important  point  is  noted.  Simple  scaling 
estimates  of  the  laminar  acoustic  boundary  layer  height  (the  Stokesian  thickness) 
5a  =  V(u/f)  indicate  that  5a  is  well  above  the  gas-phase  flame  height  except  for  very 
high  frequencies  (f^  lO'^Hz) .  Consequently,  minimal  acoustic  velocity  interaction 
with  combustion  would  be  expected  at  low  to  intermediate  frequencies,  as 
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indicated,  for  example,  by  the  more  detailed  analysis  of  Ref  6.  Although  this  does 
not  preclude  "velocity  coupled"  instability  resulting  from  laminar  interactions 
(whether  linear  or  nonlinear),  it  does  suggest  that  other  mechanisms  of 
interaction  shoiild  be  explored. 

It  is  well  established  that  for  piston  driven  closed  duct  flows  (viz.,  in  which 
there  is  negligible  mean  axial  flow),  longitudinal  acoustic  waves  of  a  few  percent 
relative  pressure  amplitude  can  induce  turbulence  within  the  Stokes  layer.  The 
literature  review  and  independent  data  of  Merkli  and  Thomann^^^  for  example, 
confirm  a  critical  amplitude  which  varies  as  Vffor  this  type  of  simple  acoustic 
motion.  Thus,  for  fixed  amplitude,  lower  frequency  disturbances  are  more  likely 
to  produce  turbtdence. 

It  is  known  from  studies  of  quasi-steady  propellant/flowfield  interactions 
that  even  small  levels  of  turbulence  within  the  combustion  zone  can  appreciably 
enhance  propellant  combustion  rates.  Whether  acoustically  induced 
turbularization  can  occur  in  actual  rocket  chambers  depends  upon  several 
complex  effects  including  the  presence  of  large  injection  rates,  surface 
Tf  ughness,  and  of  course,  the  specific  tjqje  of  combustion  process.  In  one 
experiment! 3,  values  of  the  threshold  velocity  for  propellant  response  ("acoustic 
erosivity")  were  lowered  by  a  factor  of  2  relative  to  steady  state  conditions. 
Acoustic  turbxilarization  of  the  propellant  combustion  zone  was  also  suspected  by 
Medvedev  and  Revyagin^^,  based  on  their  T-bumer  data  and  the  data  of  Price!^, 

It  is  known  from  studies  of  quasi-steady  solid  propellant/flowfield 
interactions  that  even  small  levels  of  turbulence  within  the  combustion  zone  can 
appreciably  enhance  propellant  combustion  rates,  predominantly  by  increasing 
heat  transfer  from  the  flame  zone  back  to  the  propellant  surface.  Whether 
acoustically  induced  turbularization  can  occur  in  actual  rocket  chambers 
depends  upon  several  complex  effects  including  the  presence  of  large  injection 
rates,  surface  roughness,  and  of  course,  the  specific  type  of  combustion  process. 
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Solid  rocket  motors  are  susceptible  to  several  types  of  flowfield  related 
combustion  anomalies.  Steady-state  erosion  and  acoustic  coupling  are  examples 
of  combustion  enhancements  which  can  lead  to  anomalous  burning  or,  in  some 
instances,  catastrophic  failure  of  the  motor  due  to  increased  chamber  pressures. 
The  steady-state  erosion  problem  has  been  addressed  in  past  studies  and  is 
generally  believed  to  be  a  consequence  of  the  turbulent  mean  flowfield  enhancing 
the  combustion  and  heat  transfer  processes  by  increasing  the  effective  thermal 
diffusivitylS.  Steady-state  erosive  burning  is  therefore  most  likely  to  occur  in 
rocket  motors  having  large  length  to  diameter  ratios  in  which  the  flowfields  attain 
large  velocities.  It  has  also  been  known  for  some  time  that  erosive  combustion 
can  occur  without  the  presence  of  an  appreciable  mean  flowfield,  specifically,  in 
acoustic  environments  found  in  the  center-vented  T-burner  experiments  of 
Cnimp  and  Pricers, 17  and  Medvedev  and  Revyaginl^. 

1 .2  RESEARCH  OBJECTIVES 

The  objective  of  the  present  analysis  is  to  consider,  under  several 
simplifying  assumptions,  the  general  trends  of  acoustic  boundary-layer 
turbularization  on  a  permeable  surface  in  the  presence  of  large  injection  rates. 
Though  pertinent  to  the  issues  discussed  previously,  the  nonreactive  problem  is 
more  fundamental  and  of  potentially  broader  interest.  The  analysis  is  then 
extended  to  include  the  unsteady  heat  conduction  and  combustion  of  the 
condensed  phase.  General  trends  of  reactive  acoustic  boundary-layer 
turbularization  on  a  propellant  surface,  at  the  pressure  node,  under  nominal 
homogenous  propellant  conditions  are  investigated  together  with  the  resulting 
response. 
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The  combined  response,  due  to  both  pressure  and  turbulence,  of  the  solid 
propellant  combustion  is  also  investigated.  Under  several  simplifying 
assumptions,  the  general  trends  of  the  combined  response  of  the  solid  propellant 
at  various  positions  within  a  center  vented  T-bumer,  or  low  Mach  number  region 
of  a  solid  rocket,  are  determined.  A  comparison  of  acoustic  and  steady-state 
erosion  is  also  performed. 
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SECTION  2:  ANALYSIS 


2.0  DECOMPOSITION  OF  DEPENDENT  VARIABLES 

The  motion  of  a  perfect  gas  described  by  the  Navier-Stokes  equations  is 
considered.  An  arbitrary  dependent  variable,  g(xi,t),  is  decomposed  according  to 
the  notation 

g  (Xi  ,  t)  ='g(Xi  ,  t)  +  g’(Xi  ,  t) 

with 

g  =  <g(Xi  ,  t)>  +  g"(Xi  ,  t). 

In  these  equations,  g'  is  the  turbulent  fluctuation,  g  is  the  ensemble  average,  g"  is 
the  acoustic  (deterministic)  component  and  <g>  is  the  long-time  mean.  The 
remainder  of  this  section  is  divided  into  two  parts,  the  first  considering  order-of- 
magnitude  analysis  of  the  problem  and  derivation  of  an  approximate  relation  for 
critical  amplitude  of  transition  from  the  turbulence  model  equations.  The  second 
part  discusses  the  equations  and  method  used  for  a  more  comprehensive 
numerical  solution  of  the  acoustic  boundary-layer  in  the  presence  of 
transpiration. 

U _ ORDER-OF-MAGNITUDE  ANALYSIS  18 

Neglecting  third  order  correlations  and  the  Eixial  derivatives  of  molecular 
and  turbulent  stresses,  the  <  >  average  of  the  axial  momentum  equation  is 

d 

+  —  (<p><u><u>+<p><u"u">+<u><p"u">+<u><p"u">) 

9x 

d 

—  (<p><u><v>-Kp><uV'>+<u><p"v">+<v><p"u">) 
dy 
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(1) 


9<p>  ^  d  (<n>  3  <u>)  ^  9<t> 


3x  9y  3y 


dy 


where  the  mean  turbulent  shear  stress  is  <t>  =  -<p><u'v'>  (neglecting  turbulent 

density  correlations).  The  equation  for  the  acoustic  component  (u")  is,  to  first 

order:  3u"  ,  „  „  9 ,  . 

<p> - +  (<p>u  +  <u>p  )  +  u  -^<p><u>) 

9t  9x 


3u"  .  ,..d<u> 

+  <p><u> - +  (<p>v  +  <v>p  )— — 

9x  9y 

,.9,  ,  9u"  9p’'  9,  9u"  9x'’ 

+  u  ^<p><v>)  +  <p><v> - = - + - (<|I> - )  + - 

9y  9y  9x  9y  ^  ^ 


(2) 


where  x"  =  -<p>u'v'. 

In  equation  (1),  the  second-order  correlations  arising  from  the  unsteady 
motion  can  produce  virtual  stresses  and  convective  fluxes  analogous  to  those 
produced  by  the  turbulent  motion  alone.  These  correlations  induce  the  classical 
phenomenon  of  "acoustic  streaming",l®«20  wherein  a  component  of  the  mean  flow 
is  driven  by  acoustic  motion.  Of  potentially  more  importance  to  propellant 
combustion  are  the  analogous  second-order  terms  which  appear  in  the  equation 
for  the  transport  of  thermal  energy  in  the  propellant  flame  zone.  While  these 
thermal  transport  effects  will  not  be  analyzed  in  this  study,  examples  of  their 
importance  to  effective  steady-state  heat  transfer  may  be  found  in  the  literature 
(see,  for  example.  Ref.  21).  The  effects  of  steady-state  enhancements  to  heat 
transfer  from  acoustic  streaming  effects  have  not  been  considered  with  respect  to 
the  reactive  environment  in  rocket  engines,  although  it  is  possible  that  these 
effects  could  also  produce  a  coupling  mechanism  for  instability. 

The  remainder  of  this  analysis  will  be  concerned  with  equation  (2)  and  the 
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relative  importance  of  various  terms  within  the  acoustic  boundary  layer.  To 
further  simplify  this  equation,  it  will  be  assumed  that  the  density  in  the  region  of 
interest  is  approximately  constant  and  that  the  v"  produced  by  propellant 
response  is  negligible.  These  assumptions  would  necessarily  be  removed  in 
considering  a  comprehensive  propellant  velocity  response  analysis. 

The  magnitudes  of  the  terms  in  Equation  (2)  are  estimated  by  normalizing 
the  velocities  with  Sq,  the  chamber  length  with  L,  the  chamber  radius  with  5  and 
utilizing  the  continuity  equation  and  the  assumption  of  a  standing  first 
longitudinal  acoustic  mode.  It  is  also  assumed  that  the  turbulence  shear  stress 
may  be  scaled  by  t"  =  <p>  u"^,  where  I  «  1.  The  respective  orders  of  the  terms 

in  Equation  (2)  are  then 

{ 1,  <M>,  <M>,  <M>,  (6a/8)<M>,  <M>,  (ao/f5a)<M8>,  1,  v/fsL  Cxao/(f5a,turb)  M"}. 

The  order  imity  terms  (1 )  and  (8)  recover  the  inviscid  acoustic  mode  solution 
in  the  central  region  of  the  duct.  A  considerable  simplification  results  if  <M>  is 
small,  since  terms  (2)  -  (6)  may  be  dropped  without  fundamentally  altering 
boundary  layer  characteristics  retained  in  terms  (7),  (9),  and  (10).  (It  is  important 
to  note  that  the  maximum  response  to  the  surface  sensors  in  Ref  4  occurs  in  the 
head  end  region,  where  <M>  is  ^  0.1.)  Neglecting  for  the  moment  the  transverse 
convection  and  turbulence  terms,  it  is  seen  that  for  the  viscous  term  to  be  of  order 
unity  implies  the  Stokes  estimate, 

Sa  =  (^)''2  (3) 

and  is  the  smallest  possible  laminar  boundary-layer  height  of  the  problem.  The 
effects  of  injection  on  the  acoustic  boundary  layer  [provided  largely  by  term  (7)] 
can  be  substantial  and  even  dominate.  The  ratio  between  the  convection  and 
viscous  terms  in  Equation  (3)  is  of  the  order  10  under  conditions  of  interest. 


8 


The  additional  shear  stress  provided  by  term  (10)  in  Equation  (4),  as  well  as 
the  analogous  additional  heat  flux  appearing  in  the  energy  equation,  result  from 
possible  turbulent  motions.  Further,  it  is  noted  that  term  (10)  is  the  only  retained 
term  which  depends  on  the  amplitude  of  the  acoustic  motion.  (The  acoustic 
amplitude,  M",  is  a  linear  scaling  parameter  of  all  other  terms  but  occurs 
quadraticeJly  in  the  turbulent  shear  term),  thus  posing  the  question  of  a  possible 
relation  to  nonlinear  stability  phenomena. 

As  noted  in  Ref.  12,  various  approaches  have  been  employed  to  obtain  a 
stability  or  transition  criterion  for  the  acoustic  boundary  layer.  In  the  approach 
adopted  here,  the  second-order  modeling  equations  employed  in  Ref.  5  are 
considered.  Preliminary  operations  and  assumptions  are  summarized  as 
follows.  Equations  for  the  velocity  correlations  u'iu’j  are  contracted  to  obtain  an 
equation  for  the  turbulence  intensity  q2=u'iu‘i.-  Fluid  material  properties  are 
assumed  constant,  and  the  turbulence  (or  preturbulent  disturbance  level)  is 
assumed  small  so  that  third  and  higher  order  velocity  correlations  may  be 
neglected.  It  is  also  assumed  as  an  approximation  that  u'u'  =  auuQ^» 
and  u'v'  =  auv^l^.  where  a^u,  aw  and  a^y  ®re  constants.  Substituting  these 
relations  into  the  q2  equation,  taking  the  mean  flow  to  be  negligible,  and  applying 
boimdary  layer  assumptions  for  a  quasi-planar  condition  yields: 


dq^  „  dq^  o  \ 

p(—  +  u—  +v  —  )  +  2pq(a  - +  a  - ) 

dt  dx  dy  dy 


=  -2A^  + 

a"  dy^ 


(4) 


where  A  is  the  disturbance  macrolength  -  scale  and  A  is  a  constant  used  in 
modeling  low  turbulent  Reynolds  number  dissipation.  The  approximation  A  = 
^A^a  specified  by  Equation  (3)),  is  also  employed  to  evaluate  the 

maximal  effective  length  scale  appropriate  for  the  acoustic  boundary  layer. 
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Equation  (4),  which  is  linear  in  q2,  may  be  order-of-magnitude  scaled  in  the 
same  manner  as  Equation  (2).  The  requirement  that  p(Dq2/Dt)  =  0  (for  neutral 
stability)  may  be  imposed  if  the  Reynolds  number  u"5^/u  ^  0(1).  The  order  of 

a 

magnitude  analysis  indicates  that  the  fifth  (production)  term  on  the  left  hand  side 
and  the  first  (dissipation)  term  on  the  right  hand  side  of  Equation  (7)  are 
dominant,  resulting  in  the  stability  criterion 


M"  = 


(5) 


Utilizing  the  turbulence  modeling  constants  specified  in  the  Appendix  of  Ref  3 
(i.e.,  A  =  3.25,  C\  =  .17),  and  estimating  auv  =  --15  from  fully  developed  flat-plate 
turbulent  boundary  layer  flows,  3delds  an  estimate  of  K  =  750.  Merkli  and 
Thomann^l  cite  prior  experimental  values  of  K  ranging  from  188  to  915,  and 
obtained  the  value  K  =  501  in  their  own  experiments.  The  constant  K  in  Equation 
(5)  is  related  to  the  Merkli  and  Thomann  constant  A^.  by  K  =  A^V(27t)/2.  They 
speculated  that  the  variation  in  K  observed  in  prior  studies  could  be  caused  by 
variations  in  the  roughness  of  the  duct  surface,  and  demonstrated  that 
disturbances  caused  by  the  anemometer  probe  can  also  appreciably  affect  K. 

The  critical  Mach  number  given  by  Equation  (5)  is  shown  in  Fig.  2  together 
with  the  data  of  Merkli  and  Thomann.  The  fimctional  dependence  is  correct,  and 
even  the  very  approximate  estimates  of  physical  constants  and  empirical 
parameters  yield  quantitative  agreement  to  within  several  percent  of  the  data.  To 
the  authors'  knowledge,  the  derivation  of  an  approximate  transition  relation  from 
linearizing  this  type  of  complex  turbulence  model  is  novel. 
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The  second-order  turbulence  closure  approach  developed  by  Donaldson  et 
al.22  and  Varma  et  al.23  was  implemented  in  Ref.  5  as  a  parabolized  model  for 
calculating  statistically  stationary,  compressible  transitional  flows  in  porous 
walled  ducts  with  large  injection  rates.  The  assumptions  concerning  the  gas 
phase  are  summarized  in  Appendix  A.  The  assumptions  and  order-of- 
magnitude  analysis  of  the  previous  section  again  yield  a  parabolic  equation 
system  for  the  present  acoustic  boundary-layer  problem.  Computationally,  the 
axial  convection  terms  pu  9g/9x  (where  g  is  an  arbitrary  dependent  variable)  are 
replaced  by  p9g/9t  for  the  present  case.  With  the  exception  of  9p/9x  in  the 
momentum  equation,  all  other  axial  derivatives  are  taken  to  be  null.  The 
continuity  equation. 


9t 


+ 


-V-[r''(pv  +  p'v')]  =  0 


and  the  9p/9t  term  on  the  right  hand  side  of  the  energy  equation  are  retained, 


however,  for  subsequent  reactive  flow  investigations. 


The  parabolic  differential  equation  system  may  be  considered  in  the 


functional  form 


+  pv^  =  -i-j-(r«pg^)  + 


9r  !.« 9r 


9r 


Ge(f) 


(6) 


_  .  -  ..  I  -  .  ■■■  —  m 

where  g=  { u,h,  u'u',  vV,  w'w',  uV,  h'u’,  hV,  h'h'  }  ,  u,  v  and  w  are  the  axial, 
radial  and  circumferential  velocity  components,  0  is  here  the  index  for  planar  or 
axisymmetric  geometry,  p  the  density  and  h  =  CpT  is  the  specific  static  enthalpy. 
The  molecular  transport  coefficient  |i  represents  the  dynamic  viscosity  (p)  or  the 
thermal  conductivity  parameter,  a  =  k/Cp,  as  appropriate  for  each  equation. 
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However,  not  all  of  the  molecular  diffusion  terms  for  each  equation  in  the  system 
may  be  cast  in  the  form  shown  in  Equation  (2).  Those  that  do  not  conform  are 
implicitly  contained  within  the  complex  functions  Gg,  which  also  represent  the 
sources,  cross-coupling  and  dissipation  terms  for  the  equations.  Although  the 
more  general  (  ~  )  averaging  is  used  in  the  above  notation,  the  only  mean  velocity 
component  retained  is  v. 

The  pressure  along  the  duct,  p  =  <p>  +  p"(x,t),  is  specified  by  the  one¬ 
dimensional  standing  wave  solution 

P "  =  Hma  <p>cos(n7tx/L)cos(27ift)  (7) 

where  rima  is  the  maximum  relative  acoustic  pressure  amplitude,  and 
Uma  =  Hmaat/Y.  All  nonreactive  flow  calculations  to  be  presented  were  performed 
at  the  velocity  antinode  of  the  first  longitudinal  mode.  Material  properties  were 
those  of  air  at  standard  conditions,  or  at  elevated  pressures  where  noted. 

The  turbulence  length  scale  A,  in  this  model  is  algebraic  and  described®  by 
a  linear  variation  from  its  surface  value,  Ag  (Ag  =  0  for  smooth  walls),  to  a  plateau 
level  proportionate  (with  constant  Ca)  to  the  thickness  of  the  shear  flow.  For  the 
present  analysis  this  thickness  was  defined  as  6aa>  equal  to  the  height  above  the 
surface  where  the  boundary  asymptotes  to  99%  of  u"c.  For  injected  flows, 
however,  waves  are  convected  away  from  the  surface.  In  these  cases  the  first 
zero-crossing  of  u"  was  taken  as  the  effective  5aa-  The  final  values  of  5aa  used  in 
the  length  scale  expression  were  smoothed  by  integrating  in  time.  An  alternative 
differential  length-scale  equation  would  be  more  desirable,  assuming  that 
problems  with  the  low  turbulent  Reynolds  number  behavior  of  such  equations 
could  be  circumvented. 

Boundary  conditions  at  the  duct  centerline  (or  centerplane)  are  the 
symmetry  conditions,  which  for  the  posed  system  are 
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Due  to  the  final  form  of  the  governing  equations,  the  requirement  that  v  =  0  on  the 
centerline  (centerplane)  is  necessarily  relaxed.  At  the  duct  surface,  the  static 
enthalpy  corresponding  to  a  given  temperature  is  specified,  and  the  mean 
injection  velocity  Vg  is  prescribed.  The  no-slip  condition  implies  that  all  tangential 
velocity  components  and  their  correlations  are  zero.  Experiments  have  confirmed 
that  strong  injection  rates  through  a  porous  plate  can  produce  pseudo-turbulent 
disturbances  which  must  be  included  in  the  ensemble  average  of  the  v'v' 
boundary  condition.  The  parameter  Oy  =  vV/v^  was  introduced  in  Ref.  5  to 
accoimt  for  this  important  source  of  disturbance.  All  enthalpy  correlations  are 
null  due  to  the  prescribed  iiniform  surface  enthalpy. 


The  gas  phase  equation  for  the  ensemble  mean  enthalpy,  under  the 
assumptions  outlined  in  the  previous  sections,  is  then 

+  {pv  +  +-^^(r^p  ir7) 

dt  dr  dr 

dr  \  dr  /  dt 

+  jl  1^1  +  (A  +  B  •  ReJ^  -h  Ah^  d) 

where  Ah^  co  is  the  energy  source  term  due  to  a  single-step  gas  reaction,  and  A 
and  B  are  constant  parameters  used  in  the  turbulence  model5.24  Por  the  present 
analysis  it  has  also  been  assumed  that  the  Lewis  number  is  unity,  and  that 
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combustion  of  a  homogeneous  propellant  occurs  in  a  low  Mach  number  region. 
The  turbulent  enthalpy  correlations  above  the  flame  zone  are  necessarily 
neglected  to  avoid  contragradient  diffusion  which  can  result  in  numerical 
instability.  The  species  mass  fractions  are  expressed  in  terms  of  enthalpy,  due  to 
Zeldovich  similarity24,  even  though  the  ^/3t  term  in  the  enthalpy  equation  is 
retained.  This  assumption  is  based  on  an  order  of  magnitude  analysis  of  the 
unsteady  mean  enthalpy  equation,  valid  in  the  flame  zone,  which  yields  the 
following  condition  for  similarity: 


nf « 


Pf(vs)^(cp)^ 
R  kf 


The  above  expression  shows  that  the  relative  importance  of  the  dp/dt  term  in  the 
flame  zone  is  dependent  upon  the  pressure  ratio  of  the  oscillation  and  its 
frequency  in  addition  to  the  common  parameters  associated  with  quantifying  the 
characteristics  of  the  flame  zone.  Since  the  above  expression  is  satisfied  for  the 
example  calculations  to  be  discussed  the  species  mass  fractions  needed  to 
evaluate  the  reaction  rate  expression  are  expressible  in  terms  of  enthalpy  using 
the  similarity  approach  and  the  expression  for  O)  is  then 

(H.-H)f 


ffl  =  Bg1*'exp(-TAg/T 


Ah, 


where, 

Ahg  =  I(Y„.-Yjh: 

a 


In  the  present  analysis,  ©  =  co  (t,  p)  since  it  has  been  shown  that  the  efiects  of 
near-surface  turbulent  reaction  rate  correlations  are  small  for  homogeneous 
propellants^S. 
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The  unsteady  condensed-phase  heat  conduction  is  determined  through  a 
coupled  numerical  solution  of  the  gaseous  and  condensed  phases  employing  an 
iterative  procedure  based  on  the  gas/solid  interface  energy  balance.  The 
condensed-phase  coordinate,  r\ ,  nins  from  zero  at  the  gas/solid  interface  to  -<»  in 
depth.  The  gas-phase  coordinate,  y,  runs  from  zero  at  the  gas/solid  interface  to 
+  00.  The  assumptions  concerning  the  condensed  phase  are  summarized  for 
convenience  in  Appendix  A.  The  condensed  phase  heat  conduction  equation,  in 
the  absence  of  subsurface  reactions,  is: 


aT  .  ar  a^ 


(8) 


An  energy  balance  at  the  gas/solid  interface  yields  the  boundary  condition  at  T)  =  0: 


ar 

'^an 


(9) 


where  the  enthalpy  of  decomposition  is 

Ls  =  S  (^a)g  ~ 
a 

In  equation  (9),  the  assumption  of  large  chamber  radius  is  made  for  convenience. 
Equations  (8)  and  (9)  are  both  nonlinear  due  to  the  fact  that  the  surface  regression 
rate,  r  ,  is  a  strong  function  of  the  surface  temperature,  Tg.  The  regression  rate 
(determined  by  the  condensed-phase  sublimation  or  evaporation)  is  assumed  to  be 
governed  by  an  Arrhenius  pyrolysis  law  of  the  form: 


“  =  p8V8=p7cr  =  Ase-'rAs/Tg 


(10) 
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Finally,  the  in-depth  boundary  condition  is  simply: 


Ti  ->  -oo  ,  T  -> 


(11) 


where  is  the  constant  in-depth  propellant  temperature. 

Equations  (8)  through  (11)  completely  specify  the  imsteady  conduction  of 
heat  into  the  condensed-phase.  The  solution  is  obtained  by  performing  a  linear 
coordinate  transformation  on  equations  (8)  and  (9)  to  map  the  condensed  phase 
onto  the  computational  gas  phase  grid.  Equation  (8)  and  the  gas-phase  equations 
(6)  are  then  solved  numerically,  using  the  Crank-Nicolson  finite-difference 
method,  and  equation  (9)  is  used  to  iteratively  solve  for  the  interface  temperature. 
The  linear  transformation  from  Ti-space  to  y-space  is: 


y  = 


(12) 


where  Cp  /  (k(Ts)  c^)  is  a  scaling  parameter  used  to  approximately  match 

the  depth  of  the  combustion  wave  in  the  condensed-phase  to  the  flame  height  in 
the  gas  phase.  This  transformation  enables  the  condensed-phase  solution  to  take 
full  advantage  of  the  dynamically  adaptive  gas  phase  grid  and  ensures  consistent 
numerical  treatment.  Applying  the  transformation  to  equations  (8),  (9),  and  (11) 
produces  the  following  system  of  equations: 


3T  r  dT  ^  g^i  a^T 


(13a) 


(13b) 


y->4-oo 


T^T„ 


(13c) 
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The  gas/solid  interface  energy  balance,  equation  (13b),  and  the  p3rrolysis 
relation,  equation  (10)  provide  the  coupling  between  the  condensed  and  gas  phase 
solutions.  The  coupled  solution  is  obtained  by  an  iteration  process  (here  referred 
to  as  a  substep)  at  the  new  time  step  as  follows: 

1)  Condensed  and  gas-phase  solutions  are  obtained  using  the  values  of 
the  gas/solid  interface  temperature  and  mean  injection  velocity 
specified  at  the  previous  temporal  step  or  substep 

2)  An  updated  interface  mass  flux  is  specified  using  equation  (13b)  is 
obtained  from: 

(™s)ln-i-l  ~(P7i^)ln+l 


where  n  and  n+1  are  the  n^h  and  n+1  substeps  respectively 


3)  equation  (10)  is  inverted  to  obtain  the  corresponding  surface 
temperature 


and  the  updated  gas  phase  injection  velocity  is  determined  from 
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the  condensed  and  gas  phases  are  solved  once  again,  using  the 
updated  values,  and  the  process  is  repeated  until  convergence  is 
achieved  (typically  two  or  three  iterations). 


SECTION  3:  RESULTS  AND  DISCUSSION 


3.0  STABILITY  OF  ISOTHERMAL  ACOUSTIC  BOUNDARY  LAYERS 

When  performing  isothermal  calculations,  the  initial  profile  for  the  axial 
velocity  component  was  taken  to  be  the  analytical  solution  (at  t  =  0)  of  Sexl  and 
Uchida  (in  Schlichting25)  for  oscillatory  flow  in  a  cylindrical  duct.  Initial  profiles 
of  the  normal  Re3molds  stresses  were  taken  to  be  isotropic,  and  proportionate  to 
the  square  of  the  initial  velocity  profile.  The  peak  initial  disturbance  level,  Qniax^t 
=  0,y)/ujna»  was  assumed  to  be  0.01  for  all  calculations  reported  here. 

The  implicit  numerical  procedure  and  adaptive  grid  are  described  in  prior 
studies^.  Except  where  noted,  approximately  50  time  steps  per  period  were 
employed.  The  number  of  spatial  nodes  within  the  boundary  layer  varied  from 
about  40  for  laminar  flows  to  about  75  for  turbulent  cases. 

Normalized  profiles  of  the  axial  velocity  in  the  laminar  acoustic  boundary 
layer  (f  =  1000  Hz,  6a  =  180  ^^n)  are  shown  for  various  times  within  a  period  in  Fig. 
3.  The  results  display  the  classical  phase  shift  and  Richardson  annular  effects 
due  to  viscosity.  These  calculations  were  performed  with  100  time  steps  per  period 
and  produce  a  relatively  small  computational  error,  principally  outside  the 
boundary  layer  and  near  the  maxima  in  the  acoustic  velocity.  Since  a  large 
number  of  calculations  were  needed  for  the  transition  studies,  the  number  of  time 
steps  per  period  was  halved  for  those  calculations,  and  the  error  in  the  calculated 
acoustic  velocity  was  then  approximately  doubled. 

Figure  4  shows  the  profiles  for  a  fully  (cyclically)  turbulent  boimdary  layer 
at  100  Hz  with  rTma  =  -125.  There  is  a  noticeable  diffusion  of  the  Richardson  effect 
and  a  pronounced  increase  in  the  velocity  gradient  near  the  surface.  However,  a 
protracted  law-of-the-wall  (logarithmic)  region  is  not  observed. 

Figure  5  shows  normalized  turbulence  intensity  profiles,  q/uma* 


bovindary  layer  for  the  same  conditions  stated  for  Fig.  4.  The  calculated  peaks  in 
turbulence  level  are  comparable  to  those  calculated  for  steady  state  flows  (=  1 5%), 
To  within  computational  accuracy,  these  peaks  are  symmetric  with  respect  to 
±7tfl.  A  slight  phase  lag,  on  the  order  of  5-10  degrees,  exists  between  the  intensity 
maxima  and  the  rectified  acoustic  velocity  lu"cl.  This  lag  increases  at  higher 
frequencies.  Also  note  that  the  intensity  decreases  substantially,  but  is  not 
predicted  to  vanish  at  the  zero-crossings  of  u'^. 

Figure  6  shows  turbxilent  velocity  profiles  for  a  Stokes  layer  with  injection  at 
f  =  100  Hz,  n  =  0.080,  <v>  =  1  m/s.  A  pronounced  Richardson  effect  is  evident, 
surpassing  even  the  laminar  noninjected  results.  The  estimated  value  of  the 
effective  boundary  layer  thickness  for  this  case  is  5.5  mm  (approximately  y/5a  =  14 
in  the  figure).  The  convective  wave-train  behavior  shown  is  alsv>  evident  in  the 
laminar  injected  Stokes  layer  calculations  of  Ref.  7. 

Figure  7  shows  the  corresponding  normalized  turbulence  intensity  profiles. 
Due,  in  part,  to  the  enhanced  Richardson  effect,  very  large  intensity  levels  (=  30%) 
are  predicted.  As  turbulence  is  produced  in  the  inner  layer  at  a  later  time  in  the 
cycle,  the  velocity  wave  formed  from  an  earlier  cycle  produces  a  harmonic  wave  in 
the  intensity  profile  (ft  =  5.125).  A  long  tail  is  predicted  due  to  large  decay  times 
and  the  fact  that  turbulence  is  beginning  to  accumulate  at  the  computational 
centerline  (8  =  5  cm)  where  the  symmetry  condition  is  employed.  In  an  actual 
duct  flow  environment,  the  mean  axial  flow  would  convect  these  tails  downstream 
as  they  interact  with  mean  flow  induced  turbulence. 

Calculated  results  for  the  critical  acoustic  Mach  number  at  transition, 
M"cr  =  ricr/Y*  ^^re  shown  in  Fig.  2  together  with  the  approximate  relation  and  data 
for  non-injected  Stokes  layers.  The  functional  dependence  is  similar,  but  the 
proximity  of  the  calculations  to  the  K  =  750  line  is  coincidental.  Note,  for  example, 
that  although  not  shown  here,  the  numerical  results  display  a  dependence  on  the 
initial  amplitude  of  turbulence  assumed.  "Transition"  in  the  calculated  results 
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was  determined  by  monitoring  the  growth  or  decay  of  initial  turbulence  from  the 
initial  value  over  several  cycles.  For  the  cases  with  im’ection,  the  growth  rates 
were  not  as  strongly  dependent  on  as  were  those  for  the  non-injected  cases. 

Critical  Mach  numbers  for  transition  as  a  function  of  the  normalized  mean 
injection  velocity,  <v>/(fv)l''2,  are  shown  in  Fig.  8  for  three  frequencies,  f  =  100,  300 
and  1000  Hz.  A  pronounced  decrease  in  M"cr  is  predicted  with  increasing 
injection  velocity,  so  that  all  curves  exhibit  a  minimum  that  is  frequency 
dependent.  The  minima  are  almost  a  factor  of  three  below  the  M"cr  for 
noninjected  Stokes  layers.  As  an  example,  the  100  Hz  non-injected  Stokes  layer  is 
calculated  to  have  M'^j.  =  .083.  The  minimum  value  of  M"cr  =  025  occurs  at  <v>  = 
0.25  m/s,  while  M"cj.  =  .048  at  <v>  =  1.0  m/s  (the  termination  of  the  line  at  this 
frequency). 

Even  more  informative  are  the  same  results  for  transition  normalized  not 
as  a  critical  acoustic  Mach  number,  but  as  a  critical  acoustic  Reynolds  number, 
UniaAfu)^'^^.  In  this  case,  all  three  curves  for  frequency  in  Fig.  8  collapse  onto  a 
single  line  (to  within  computational  precision),  shown  as  the  upper  line  in  Fig.  9. 
This  figure  also  shows  results  for  finite  levels  of  the  surface  disturbance 
parameter  Oy  =  0.035  and  0.07  for  two  different  transition  criteria  to  be  discussed. 
The  accompanying  value  of  Ag  assumed  in  the  calculations  was  3x1 0’^  m.  These 
values  are  viewed  as  realistic  and  perhaps  conservative.  It  is  not  surprising  that 
this  finite  disturbance  effect,  which  continuously  "feeds"  the  injected  layer,  is 
predicted  to  3neld  further  reduction  in  critical  acoustic  Reynolds  number, 
analogously  to  quasi-steady  flows.^  As  a  specific  example,  the  minimum  critical 
Mach  number  for  the  100  Hz  case  at  <v>  =  .20  m/s  is  decreased  to  approximately 
one  percent  for  Oy  =  0.035. 

As  mentioned  earlier,  it  was  more  difficult  to  determine  the  value  of  M"ci. 
for  cases  with  injection  and  finite  levels  of  Oy  present.  The  first  criterion  employed 
consisted  of  monitoring  the  maximum  turbulence  level  qmax^^ma  occurring  at 


21 


any  height  within  the  boundary  layer.  If  this  maximum  attained  a  value  of  0.07 
or  greater  within  3  cycles,  the  flow  was  considered  turbulent.  Using  this 
criterion,  the  Oy  =  0  lines  were  computed  along  with  the  altemating-dashed  lines 
for  Oy  =  0.035  and  0.07.  While  this  criterion  was  satisfactory  for  Oy  =  0  cases,  it 
proved  to  be  unsatisfactory  for  cases  with  finite  levels  of  Oy. 

To  illustrate  this  effect.  Fig.  10  shows  the  variation  of  turbulence  intensity 
vs  normalized  time  at  y  =  6a,  for  three  values  of  U,  <v>  =  0.25  m/s  and  Oy  =  0.  It  is 
evident  from  the  cvirves  in  Figs.  10a  and  10b  that  a  very  small  increase  in  IT  is  able 
to  cause  the  flow  to  transition  to  a  cyclically  turbulent  condition. 

When  finite  levels  of  Cy  are  included,  the  determination  of  Her  becomes 
more  complicated.  Figures  11a,  b,  and  c  are  results  for  <v>  =  0.25  m/s  and  Oy  = 
0.035.  Figure  11a  is  a  case  for  which  the  turbulence  bursts  to  a  q/Uj^j^  level  of 
approximately  0.24  but  then  dissipates  over  10  cycles  imtil  only  the  low  level 
turbulence  due  to  Oy  is  still  present.  Figures  11b  and  c  illustrate  the  increase  in  fl 
required  to  achieve  the  cyclic  behavior  which  we  considered  to  be  characteristic  of 
a  turbulent  case.  This  criterion  for  transition  was  used  to  calculate  the  dashed 
curves  for  and  0.07  in  Fig.  9. 

To  recapitulate  the  results  of  the  two  transition  criteria,  the  alternating 
dashed  lines  in  Fig.  9  qualitatively  represent  a  "burst"  or  transient  criterion  for 
transition,  while  the  dashed  lines  represent  the  particular  (long-time)  criterion. 


For  the  latter  case,  regions  to  the  right  of  the  nearly  vertical  line  segments  are 
indicated  to  be  turbulent.  The  difference  between  the  two  criteria  is  in  part 


explained  by  Equation  (4),  which  must  be  augmented  by  the  term  p  <v>  when 
is  finite.  This  can  result  in  three  different  turbulence  time  scales  being 


competitive  under  some  conditions. 
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The  therinochemical  properties  and  reference  operating  conditions  of  the 
baseline  solid  propellant,  presented  in  Appendix  B,  were  chosen  to  represent  the 
nominal  properties  of  homogeneous  propellants  and  were  also  used  in  the  erosive- 
burning  study  of  Ref  26.  Calculations  presented  were  performed  using  the 
baseline  propellant,  which  has  an  exothermic  surface  reaction,  and  an  identical 
propellant  with  an  endothermic  surface  reaction,  L°=837  J/kg.  To  reproduce  the 
reference  operating  conditions,  the  pre-exponential  constant  in  the  gas  phase 
reaction  rate  was  calculated  to  be  8.2x109  kgi  *^m‘**“^s'lKBg  for  the  endothermic 
propellant.  Steady-state  calculations  performed  over  a  broad  range  of  pressures 
yielded  a  normal  rate  of  regression  r  =  f  (p/p*)"  'vit^  r  and  n  as  indicated  in  the 
Table.  Transition  calculations  were  performed  using  50  time  steps  per  cycle  and 
two  iterations  per  step  while  propellant  response  calculations  were  performed 
using  100  steps  per  cycle  and  three  iterations  per  step.  Unless  stated  otherwise, 
all  propellant  responses  were  obtained  at  an  acoustic  frequency  of  1000  Hz  and  five 
percent  acoustic  pressure  ratio,  Hma  =  0.05.  The  disturbance  level  of  the  gas 
transpiring  from  the  propellant  surface,  Oy,  was  assumed  to  be  zero  in  all 
calculations. 

Shown  in  Fig.  12  are  the  normalized  steady-state  temperature  profiles 
obtained  in  the  gas  and  condensed  phases  for  the  endothermic  and  exothermic 
propellant  formulations  at  reference  conditions,  normalized  by  the  approximate 
theoretical  gas  phase  flame  height,  5f,  (note  that  the  approximate  5f  scaling  used 
here  is  independent  of  L°).  The  results  indicate  that  the  exothermic  formulation 
provides  an  anticipated  increased  gas  phase  flame  height,  increasing  the 
possibility  of  acoustically  induced  turbulence  interacting  with  the  flame  zone. 
The  coincident  condensed  phase  profiles  are  shown  plotted  in  the  gas  phase 


computational  coordinates  and  are  therefore  scaled  by  ^^  =  4.0  in  all  calculations. 
Thus,  if  the  scaling  were  not  employed,  the  condensed  phase  combustion  wave 
would  be  approximately  four  times  the  length  of  the  flame  zone  and  the  condensed 
phase  solution  would  not  efficiently  use  the  dense  grid  (generally  35  points) 
provided  by  the  dynamically  adaptive  grid  within  the  combustion  zone. 

Figure  13  presents  the  normalized  temperature  and  turbulent  velocity 
profiles  vs.  distance  from  the  burning  surface.  The  abscissa  is  normalized  by 
5a  =  100  pm.  The  profiles  are  given  for  various  times  within  the  acoustic  cycle  for 
the  baseline  propellant  with  an  equivalent  roughness  of  kg  =  1 50  pm.  The  abscissa 
is  presented  as  a  logarithmic  scale  to  reveal  the  details  in  the  near  surface  flame 
region.  The  results  obtained  show  the  phase  shift  and  Richardson  annular  effect 
due  to  viscosity  and  injection  that  were  obtained  in  the  previous  section  for 
isothermal  Stokes  layers.  The  acoustic  velocity  is  negligible  in  the  region  of  the 
flame  zone  at  this  intermediate  frequency.  However,  Fig  14  shows  the 
corresponding  normalized  temperature  profiles  and  turbulence  intensity  profiles. 
The  peaks  in  turbulence  intensity  within  the  cycle  correspond  approximately  to 
the  large  velocity  gradients  produced  by  the  Richardson  annular  effect.  The 
results  indicate  that  the  turbulence  intensity  obtains  magnitudes  of  one  to  two 
percent  within  the  combustion  region  at  certain  points  within  the  acoustic  cycle, 
and  results  in  propellant  combustion  responding  to  enhanced  heating  rates  due  to 
turbulent  diffusion. 

In  Section  3.0,  it  was  found  that  for  a  smooth  surface  (kg  =  0)  without 
surface  injection  disturbances  (Oy  =  0),  an  acoustic  transition  locus  could  be 
expressed  in  terms  of  axial  acoustic  velocity  Reynolds  number  as  a  function  of  an 
injection  velocity  acoustic  Reynolds  number.  Figure  15  shows  the  previously 
reported  results.  The  critical  amplitude  locus  is  quite  sensitive  to  the  method 
used  to  introduce  "seed"  turbulence  into  the  boundary  layer  and  to  other 
conditions  such  as  the  relative  height  within  the  boundary  layer  at  which 


amplification  is  determined.  Figure  15  shows  additional  calculations  for  an 
isothermal  boimdary-layer  which  utilize  an  alternative  means  of  introducing  the 
initial  disturbance.  (The  difference  between  the  two  methods  is  that  the  earlier 
results  introduced  a  1%  peak  disturbance  level  within  the  Stokes  layer  thickness 
5a,  while  the  present  results  introduce  the  same  1%  disturbance  level  over  the 
much  larger  actual  injected  boundary  layer  thickness,  5aa)- 

Using  this  new  and  perhaps  more  realistic  method,  transition  calculations 
were  also  performed  for  the  baseline  propellant  at  frequencies  of  100  and  1000  Hz 
over  pressures  ranging  from  10  to  100  atmospheres.  The  transition  points 
obtained  for  propellant  combustion  are  in  good  agreement  with  the  isothermal 
calculations.  Thus,  the  appropriate  scaling  involves  flame  or  edge  quantities 
only,  and  indicates  that  the  acoustic  flowfield  (for  a  smooth  surface)  can  be 
approximated  as  an  isothermal  flow  at  the  steady-state  flame  temperature  in 
determining  transition. 

Figure  16  shows  the  calculations  presented  in  Fig  15,  for  the  endothermic 
propellant  along  with  calculations  for  the  exothermic  propellant  at  an  acoustic 
frequency  of  100  Hz,  plotted  as  a  function  of  the  critical  acoustic  pressure  ratio, 
Her,  vs.  the  mean  pressure  normalized  by  p*.  The  critical  pressure  ratio  is  foimd 
to  be  insensitive  to  the  heat  released  from  the  surface  reaction,  although  the 
actual  response  of  the  propellant  is  dependent  on  L®.  This  results  from  the  fact 
that,  for  a  smooth  wall,  transition  is  initiated  above  the  thin  combustion  zone. 
The  critical  pressure  ratio  is  found  to  decrease  with  increasing  mean  pressure  in 
all  cases,  indicating  that  a  propellant  which  may  be  stable  with  respect  to  a 
specific  acoustic  mode  at  a  known  mean  pressure  may  encounter  acoustic 
transition  at  that  mode  when  operating  at  an  elevated  mean  pressure.  This 
phenomena  may  explain  stability  problems  encountered  in  the  development  of 
solid  rocket  motors  with  unconventionally  high  chamber  pressures. 

Figure  1 7  shows  the  response  of  the  baseline  endothermic  propellant  to  the 
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presence  of  the  turbulent  acoustic  boundary-layer  for  various  surface 
roughnesses.  The  normalized  propellant  mass  flux  (burn  rate)  is  plotted  vs. 
normalized  time,  x.  A  diagram  of  the  variation  of  the  acoustic  edge  velocity  and 
pressure  is  included  at  the  top  of  the  figure.  The  pressure  variation  is  correct  for 
an  axial  location  in  the  head-end  region  of  a  chamber.  In  the  afl-end  region,  the 
pressure  would  be  shifted  180  degrees  in  phase  with  respect  to  the  pressure  curve 
shown.  The  case  of  a  smooth  wall,  kg=  0  m,  required  a  15  percent  pressure  ratio 
to  obtain  the  response  shown.  The  cases  in  which  a  finite  value  of  surface 
roughness  was  prescribed  were  calculated  using  a  five  percent  pressure  ratio. 
The  smooth  wall  produced  only  a  slight  response  in  the  propellant  burn  rate  that 
is  rectified  with  respect  to  the  velocity  field  and  which  is  negligible  over  a  large 
portion  of  each  half  cycle.  The  calculations  in  which  finite  levels  of  surface 
roughness  are  included  reveal  several  interesting  features.  All  cases  exhibit  the 
rectified  response  along  with  a  time-mean  augmentation  ("D-C  shift")  of  the 
propellant  burn  rate,  about  eight  percent  for  kg  =  150pm,  five  percent  for  125pm, 
and  two  percent  for  the  100pm  surface  roughness.  The  relative  phase  of  the 
response,  with  respect  to  the  acoustic  velocity  and  pressure  field,  varies  with 
surface  roughness  and  may  provide  a  mechanism  by  which  the  acoustical  energy 
can  be  reinforced  or  damped  according  to  the  Rayleigh  criterion. 

Figure  18  compares  the  effect  of  surface  decomposition  energy  on  the 
propellant  response.  The  calculations  were  performed  for  the  endothermic  and 
exothermic  propellant  formulations  with  kg  =  150pm,  and  n  =  0.05,  in  both  cases. 
The  exothermic  propellant  has  a  mean  shift  in  burning  rate  that  is  approximately 
three  times  as  large  as  that  occurring  in  the  endothermic  case. 

The  propellant's  response  to  the  acoustic  boundary  layer  is  further 
investigated  using  a  conventional  Fast  Fourier  Transform  (FFT)  analysis  over 
several  cycles.  The  cases  examined  are  identical  to  the  calculations  with  finite 
roughness  presented  in  Fig  17  and  discussed  above.  Figures  19  and  20  present  the 
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results  of  the  analysis  in  terms  of  the  relative  amplitude  and  frequency  content  of 
the  response  and  the  relative  phase  of  the  responses  at  each  frequency, 
respectively.  The  large  dashed  curve  is  the  driving  frequency  of  the  external 
acoustic  velocity,  the  small  dashed  curve  is  the  response  of  the  propellant  with  a 
smooth  wall  (the  response  has  been  graphically  shifted  +300  Hz  to  facilitate  the 
comparison),  and  the  solid  curve  is  the  response  with  kg  =  150  ^m.  The  frequency 
content  and  relative  amplitudes  of  the  response  with  kg  =  150  pm  is  representative 
of  the  cases  with  finite  surface  roughness.  The  results  presented  in  Fig  19  have, 
in  each  case,  been  normalized  by  the  largest  amplitude  obtained  for  that 
calculation  after  removal  of  the  zero-frequency  component.  It  can  be  seen  that  in 
both  calculations  the  primary  response  occurs  at  the  second  harmonic  with 
significant  harmonic  content  occurring  at  the  lower,  even  harmonics.  Only  even 
harmonics  are  obtained  for  this  case  since  axial  mean  flow  is  neglected  and 
calculations  are  performed  at  a  pressure  node.  The  relative  phase  of  the 
propellant  response,  presented  in  Fig  20,  is  seen  to  vary  substantially  from 
harmonic  to  harmonic.  It  can  further  be  seen  that  at  the  second  harmonic,  2000 
Hz,  the  phase  lag  of  the  response  increases  with  increasing  surface  roughness 
(excluding  the  smooth  wall  case  due  to  the  large  pressure  ratio  required  to  obtain 
the  response). 

3.2  COMBINED  RESPONSE  AND  ACOUSTIC  EROSION  AT  VARIOUS 
AXIAL  LOCATIONS 

The  thermophysical  properties  and  reference  operating  conditions  of  the 
solid  propellant  employed  in  the  previous  section  were  also  used  to  investigate  the 
combined  response  and  acoustic  erosion  at  various  axial  locations.  All 
calculations  presented  were  performed  using  100  time  steps  per  cycle  and  three 
iterations  per  time  step.  When  calculating  the  acoustic  erosion  of  the  propellant 


27 


at  various  axial  locations,  the  integration  of  the  propellant  response  was 
performed  in  the  7th  cycle  of  oscillation  after  the  introduction  of  the  seed 
turbulence  in  order  to  let  the  propellant  combustion  reach  a  cyclically  steady- 
state.  The  integration  required  to  calculate  the  acoustic  erosion  due  to  acoustic 
velocity  oscillations  alone  were  performed  in  the  3rd  cycle  of  oscillation  after  the 
seed  turbulence  was  first  introduced  since  the  flow  would  reach  the  cyclically 
steady-state  condition  in  approximately  two  cycles.  All  calculations  presented 
were  performed  at  an  acoustic  frequency  of  1000  Hz  and  five  percent  acoustic 
pressure  ratio,  Hma  =  0.05,  unless  otherwise  stated.  The  calculated  acoustic 
boundary  layer  heights  were  8aa  =  1mm  and  0.6mm  for  an  acoustic  frequency  of 
f=1000  Hz  and  2000  Hz,  respectively. 

The  condensed-phase  solution  was  calculated  on  the  gas-phase 
computational  grid  after  being  scaled  by  ^^  =  4.0  to  use  the  dense  grid  provided  by 
the  combustion  zone.  The  flame  zone,  which  has  a  height  of  approximately  lOpm, 
generally  had  35  grid  points  within  it  while  the  acoustic  boundary-layer  and 
remainder  of  the  computational  domain,  which  extends  to  a  height  of  5cm, 
contained  approximately  60  and  180  grid  points,  respectively.  The  disturbance 
level  of  the  gas  transpiring  from  the  propellant  surface,  Oy,  was  assumed  to  be 
zero  in  all  calculations.  The  acoustic  erosive  burning  is  calculated  by  integrating 
the  instantaneous  propellant  regression  rate,  r,  over  one  cycle  of  oscillation  to 
obtain  the  long  time  mean  of  the  regression  rate,  <  r  >,  which  is  then  normalized 
by  the  normal  burning  rate  of  the  propeUant,  rp. 

Shown  in  Fig.  21  are  the  steady-state  temperature  profiles  in  the  gas  and 
condensed-phases  normalized  by  the  steady-state  flame  temperature,  Tf,  as  a 
function  of  distance  from  the  burning  surface  normalized  by  the  approximate 
theoretical  laminar  flame  height,  5f.  The  condensed-phase  temperature  profile  is 
plotted  in  the  gas-phase  computational  coordinates  and  has  therefore  been  scaled 
by  Cn  =  4.0  as  it  is  in  all  calculations.  The  calculated  profiles  are  in  good 
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agreement  with  the  approximate  theoretical  flame  height  and  theoretical 
condensed-phase  thermal  boundary-layer  depth. 

Figure  22  is  a  comparison  of  the  gas-phase  temperature  perturbations  due 
to  acoustic  pressure  oscillations  (laminar  calculation)  calculated  with  and 
without  the  assumption  of  Zeldovich  similarity.  The  latter  case  being  equivalent 
to  the  equation  system  used  in  Tien’s  analysis. 2 7  The  calculations  were 
performed  at  x/L=0.375  with  nma=0.05  (which  yields  a  local  acoustic  pressure 
ratio  of  1.913  percent  (ni=0.01913)  and  are  foimd  to  be  virtually  identical.  The 
perturbations  presented  are  defined,  following  the  analysis  of  Ref.  27,  as  the 
difference  between  the  instantaneous  value  and  the  steady-state  value  of  the 
variable  normalized  by  a  reference  condition  (normally  a  value  occurring  at  the 
edge  of  the  flame  zone)  and  the  local  acoustic  pressure  ratio.  Various  profiles 
throughout  one  cycle  of  oscillation  are  plotted  as  a  function  of  normalized  distance 
from  the  burning  surface  to  a  distance  of  1000  theoretical  flame  heights.  The 
perturbations  qualitatively  agree  with  the  results  of  Ref.  27  and  veiify  the 
oscillatory  behavior  of  the  wave  train  exterior  to  the  flame  zone.  A  direct 
comparison  with  the  results  of  Ref.  27  was  not  possible  since  thermochemical 
parameters  of  a  representative  solid  propellant  were  selected  for  this  study.  This 
would  imply  that  the  (dimensional)  frequencies  required  for  a  direct  comparison 
with  the  results  of  Ref  21  be  the  order  of  10®  Hz.  At  frequencies  of  this  magnitude, 
many  of  the  assumptions  of  both  analyses  would  have  to  be  reevaluated. 
Calculations  up  to  10'^  Hz  were  performed  which  appeared  to  be  in  qualitative 
agreement  with  the  perturbation  results  of  Ref  27. 

The  curves  presented  in  Fig.  23  are  of  the  laminar  and  turbulent  gas-phase 
temperature  perturbations  occurring  at  various  points  within  one  cycle  of 
oscillation  plotted  as  a  function  of  normalized  distance  from  the  burning  surface 
corresponding  to  the  gas-phase  flame  zone  plotted  in  Fig.  21.  The  calculations 
were  performed  at  x/L=0.375,  with  ni=0.01913  and  f=1000  Hz.  The  laminar 
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perturbations  oscillate  symmetrically  about  the  steady-state  solution,  as  expected, 
since  the  pressure  oscillates  sinusoidally.  The  turbulent  perturbations  are  found 
to  all  be  positive,  accounting  for  the  increase  of  approximately  5  percent  in  the 
mean  regression  rate  of  the  condensed-phase  (DC-shifl)  at  this  location,  and  are 
substantially  larger  in  magnitude  than  the  corresponding  laminar  perturbations. 
It  is  noted  that  the  "turbulent"  perturbations  presented  are  actually  the  combined 
response  due  to  both  turbulent  flow  and  harmonic  pressure  fluctuations,  but  will 
be  referred  to  as  the  "turbulent"  perturbations  in  subsequent  text  for  brevity. 

Figure  24  compares  the  laminar  and  turbulent  gas  phase  normal  velocity 
perturbations  directly  corresponding  to  the  temperature  perturbations  presented 
in  Fig.  23.  The  laminar  perturbations  have  been  scaled  in  magnitude  by  a  factor 
of  5  for  presentation.  It  is  evident  that  the  pressure  response  is  dominated  by  the 
effects  of  turbulence  produced  by  the  acoustic  motion.  The  turbulent  normal 
velocity  perturbations  at  the  surface  and  edge  of  the  flame  zone  are  found  to 
oscillate  about  a  mean  shift  corresponding  to  the  shift  in  the  regression  rate  of  the 
condensed-phase. 

The  laminar  and  turbulent  perturbations  of  the  condensed-phase 
temperature,  directly  corresponding  to  the  curves  presented  in  Figs.  21,  23  and 
24,  and  plotted  in  the  gas-phase  computational  space  are  shown  in  Fig.  25.  The 
magnitude  of  the  laminar  perturbations  have  been  scaled  by  a  factor  of  5  for 
presentation.  The  laminar  perturbations  do  not  oscillate  about  the  steady-state 
profile,  as  the  laminar  gas-phase  perturbations  do,  but  oscillate  about  an 
asymmetrical  mean  perturbation.  This  asymmetry  arises  from  the  variation  of 
the  regression  rate  (and  therefore,  convection)  of  the  condensed-phase  throughout 
the  cycle  of  oscillation  in  combination  with  the  oscillatory  nature  of  the 
temperature  at  the  gas-solid  interface.  The  turbulent  perturbations  are  foxmd  to 
be  on  the  order  of  ten  times  larger  than  the  corresponding  laminar  perturbations. 
The  extremely  as3mmetrical  behavior  of  the  turbulent  profiles  is  due  primarily  to 
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the  increase  in  the  mean  regression  rate  of  the  condensed-phase  at  this  location, 
which  tends  to  reduce  the  depth  of  the  thermal  boundary-layer.  The  laminar  and 
turbulent  acoustic  velocity  profiles  and  the  corresponding  turbulence  intensity 
profiles  are  essentially  identical  to  the  corresponding  figures  of  the  previous 
section. 

Results  have  been  obtained  for  the  response  of  the  propellant  regression 
rate  to  conditions  which  would  be  representative  of  those  existing  in  a  center- 
vented  T-bumer,  or  in  the  head-end  region  of  a  solid  rocket  motor  where  the  mean 
axial  Mach  number  is  low,  when  either  are  experiencing  acoustic  oscillations. 
For  the  fundamental  mode,  the  acoustic  environment  varies  from  the  pure 
pressure  oscillation  at  either  end  to  the  pure  velocity  oscillations  in  the  center  of 
the  chamber  at  the  pressure  node,  and  a  combination  of  both  at  any  point  between. 
Only  results  for  the  portion  of  the  duct  from  x/L=0  to  x/L=0.5  will  be  presented 
since  in  the  absence  of  a  mean  axial  flow  the  response  will  be  cyclically  identical 
in  either  end  of  the  duct. 

The  combined  response  of  the  condensed  phase  regression  rate  at  various 
positions  along  the  chamber  are  shown  in  Fig.  26.  The  top  portion  of  the  plot 
contains  a  diagram  of  the  corresponding  acoustic  pressure  oscillation  and  the  two 
possible  external  acoustic  velocity  oscillations  which  could  exist  at  each  location 
along  the  chamber  (the  diagram  shows  only  the  relative  phases  of,  and  is  not 
intended  to  show  the  amplitude  relationship  between,  the  acoustic  velocity  and 
pressure).  The  calculations  were  performed  for  Fima  =  0.05  and  kg  =  150|im.  The 
response  of  the  condensed-phase  regression  rate  at  the  end  of  the  chamber,  x/L=0, 
is  sinusoidal  and  lags  the  pressure  by  approximately  9  degrees.  At  a  quarter 
point  in  the  duct,  x/L=0.25,  the  turbulent  acoustic  boundary-layer  produces  a 
slight  increase  in  the  mean  regression  rate  (1 .8  percent)  on  which  the  sinusoidal 
response  of  the  regression  rate  to  the  acoustic  pressure  is  superimposed.  As  the 
location  is  moved  closer  to  the  center  of  the  chamber,  nearer  the  pressure  node. 
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the  acoustically  produced  turbulence  begins  to  dominate  the  dynamics  of  the 
combustion  zone.  At  x/L=0.3125  the  sinusoidal  behavior  of  the  response  has 
nearly  vanished  as  the  response  becomes  increasingly  rectified.  At  x/L=0.4375  the 
response  is  almost  completely  rectified  due  to  the  turbulence,  yet  the  small 
acoustic  pressure  oscillation  which  exists  produces  discemable  variations  in  the 
peaks  and  troughs  of  the  response.  At  the  pressure  node,  x/L=0.5  (not  shown),  the 
response  is  completely  rectified  and  has  a  mean  shift  of  6.6  percent,  the 
maximum  in  the  chamber. 

The  response  of  the  condensed-phase  regression  rate  to  the  acoustic 
boundary  layer  is  further  investigated  using  a  conventional  Fast  Fourier 
Transform  (FFT)  analysis  on  6  cycles  of  oscillation  containing  600  points  of  data. 
The  variation  of  the  frequency  content  of  the  response  along  the  axis  of  the 
chamber  was  calculated  for  f=1000  Hz,  Tima  =  0.05,  and  two  equivalent  sand 
roughness  heights,  ks=150pm  and  ks=100)im.  Fig.  27  :.hows  the  normalized 
frequency  content  for  ks=l  50pm  at  various  locations  along  the  duct  and  includes 
the  responses  plotted  in  Fig.  26.  The  results  presented  at  each  axial  location  have, 
in  each  case,  been  normalized  by  the  largest  amplitude  obtained  in  the  entire 
chamber  after  removal  of  the  zero  frequency  component  corresponding  to  the  DC- 
shift  in  the  mean  regression  rate.  The  primary  response  from  x/L=0  to  x/L=0.25  is 
seen  to  be  that  of  the  driving  frequency  corresponding  to  the  acoustic  pressure 
oscillation.  The  lack  of  any  harmonic  content  in  this  region  is  due  to  the  nearly 
first  order  (0=1.1)  gas-phase  reaction.  In  the  region  from  x/L=0.3125  to  the  center 
of  the  chamber  the  frequency  content  is  characterized  by  a  steady  decrease  in  the 
amplitude  of  the  fundamental  frequency  and  a  rapid  rise  in  the  amplitude  of  the 
second  harmonic  corresponding  to  the  rectified  turbulent  response.  There  is  a 
similar  appearance,  and  growth  in  the  amplitudes  of,  the  higher  even  harmonics 
along  with  the  appearance  of  some  very  small  odd  harmonics.  The  largest 
response  is  seen  to  occur  at  the  center  of  the  chamber  due  to  the  acoustically 
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induced  turbulence  and  is  approximately  double  the  amplitude  of  the  largest 
pressure  response,  which  occurs  at  x/L=0.  Examining  the  frequency  content  at 
the  center  of  the  chamber  reveals  that  only  a  completely  rectified  response  with 
even  harmonics  is  produced,  as  was  shown  in  the  previous  section,  by  the 
turbulent  acoustic  velocity  oscillations  alone.  It  can  be  seen  that,  although  no  odd 
harmonics  are  produced  by  either  the  pure  acoustic  pressure  or  turbulent 
acoustic  velocity  oscillations  alone,  in  combination  they  do  produce  detectable  odd 
harmonics.  The  calculated  FFT  results  corresponding  to  the  100pm  equivalent 
sand  roughness  height  were  qualitatively  identical  to  the  results  presented  in  Fig. 
27,  differing  only  in  the  magnitude  of  the  turbulent  response. 

The  results  presented  in  Figs.  26  and  27  may  also  be  examined  to  assess 
whether  the  acoustic  oscillations  within  the  chamber  could  be  damped  or 
reinforced  by  the  propellant  response.  The  frequency  response  presented  in  Fig. 
27  indicates  that  an  acoustic  oscillation  at  the  fundamental  frequency  is  capable  of 
introducing  energy  to  acoustic  modes  at  the  higher  harmonic  frequencies.  The 
response  of  Fig.  26  indicates  that  the  energy  added  over  one  cycle  to  the  acoustic 
oscillations  would  be  negligible  towards  the  center  of  the  chamber  due  to  the 
nearly  complete  rectification  of  the  combined  response.  The  potential  for 
combustion  instability  arising  in  an  actual  rocket  chamber,  due  to  the  combined 
response,  will  depend  upon  the  extent  to  which  the  turbulent  component  is  biased, 
and  shifted  in  phase,  with  respect  to  the  local  pressure  oscillation  by  the 
interaction  of  the  developing  mean  axial  flow  with  the  acoustic  boundary-layer 
and  turbulence  fields.  The  relative  magnitude  of  the  turbulent  and  pressure 
components  of  the  combined  response  indicate  that  turbulent  velocity  coupling 
could  potentially  contribute  a  larger  portion  of  acoustic  energy  to  the  overall 
flowfied  than  that  due  to  pressure  coupling  alone. 
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The  experiments  of  Crump  and  Priced  5,1 7  using  center- vented  T-bumers  to 
investigate  the  combustion  of  solid  propellants  in  an  oscillatory  environment, 
produced  valuable  data  concerning  the  affects  of  acoustic  erosion  throughout  the 
chambers  of  such  devices.  Ref.  15  presents  data  on  the  variation  of  acoustic 
erosion  throughout  the  length  of  the  chamber  of  the  T-bumer.  The  results  all 
exhibit  a  maximum  increase  in  mean  regression  rate  at  the  pressure  node, 
corresponding  to  the  maximum  in  acoustic  oscillations,  and  a  minimum  (always 
negative)  in  mean  regression  rate  at  each  end  of  the  chamber.  They  inferred  that 
the  erosion  was  produced  by  a  negative  pressure  response  near  each  pressure 
antinode  and  a  positive  velocity  response  near  the  pressure  node.  The  Soviet 
experimenters,  Medvedev  and  Revyaginl^,  verified  the  results  of  Crtunp  and  Price 
with  T-bumer  experiments  of  their  own  and  suggested  that  the  large  increase  in 
the  mean  regression  rate  may  be  due  to  the  turbularization  of  the  Stokes  layer 
close  to  the  surface  of  the  burning  propellant.  Fig.  28  shows  the  variation  of  the 
acoustic  erosion  as  a  function  of  axial  location  within  the  chamber  for 
calculations  with  ks=100p.m,  125)im,  and  150ji^.  All  of  the  calculations  were 
performed  at  f=1000  Hz  and  lima  =  0-05.  An  equivalent  sand  roughness  height  of 
150p,m  is  found  to  produce  an  acoustic  erosion  of  approximately  1  percent  from 
very  near  the  end  of  the  chamber  to  x/L=0.3125,  where  the  erosion  begins  to  rise 
rapidly  to  a  maximum  of  6.6  percent  at  x/L=0.5.  Values  of  lOO^im  and  125^m 
equivalent  sand  roughness  heights  produce  similar  variations  of  erosion 
throughout  the  chamber,  differing  substantially  only  in  their  respective 
magnitudes.  The  acoustic  erosion  remains  zero  to  approximately  x/L=0.25  and 
then  increases  to  the  maximum  values  of  1 .6  percent  and  4  percent  at  x/L=0.6. 
The  lack  of  negative  erosion  near  the  pressure  antinodes  in  the  present 
calculations  may  be  a  result  of  the  choice  of  propellant  thermophysical 
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parameters  or  the  simplifying  assumptions  made  for  the  combustion  process. 

Fig.  29  shows  additional  experimental  results  concerning  erosive  burning 
obtained  by  Crump  emd  Price  and  Medvedev  and  Revyagin  as  presented  in  Ref.  14. 
The  mean  shift  in  the  condensed  phase  regression  rate  is  plotted  as  a  fimction  of 
the  velocity,  Ue,  for  both  acoustic  and  steady-state  erosive  conditions.  Here,  Ug 
represents  the  centerline  velocity  of  the  steady-state  experiments  and  the 
amplitude  of  the  acoustic  velocity  oscillation  in  the  oscillatory  experiments.  A 
different  propellant  was  used  by  each  set  of  experimenters,  but  the  same  trends 
were  discerned.  First,  the  acoustic  erosion  threshold  velocity  was  found  to  be 
approximately  half  of  the  threshold  velocity  for  steady-state  erosion.  Second,  the 
acoustic  erosion  was  noted  to  be  approximately  linearly  related  to  the  maximum 
amplitude  of  the  acoustic  velocity  oscillation. 

Fig.  30  shows  results  of  calculations  of  acoustic  and  steady-state  erosion  for 
the  propellant  considered  in  the  present  investigation,  assuming  two  values  of 
equivalent  sand  roughness  height,  ks=150pm  and  100pm.  The  acoustic 
calculations  were  performed  at  the  pressure  node  to  insure  that  the  response  was 
due  solely  to  the  acoustically  induced  turbulence.  The  maximum  acoustic  velocity 
amplitude  a  calculation  can  be  performed  for  is  limited  by  the  corresponding 
maximum  allowable  Reynolds  niunber,  based  on  wall  shear  stress  and  equivalent 
sand  roughness  height,  above  which  the  validity  of  the  surface  roughness  model 
is  questionable26.The  steady-state  calculations  were  performed  using  the  SPEC24 
code  at  the  same  operating  conditions,  and  for  the  same  propellant  formulation, 
as  the  acoustic  calculations.  The  present  calculations  predict  the  acoustic 
threshold  velocity  to  be  less  than  half  the  threshold  velocity  of  steady-state  erosion 
and  the  acoustic  response  is  seen  to  vary  nearly  linearly  with  velocity  after  the 
initial  response  appears.  Both  of  these  results  agree  well  with  the  experimental 
results  shown  in  Fig.  29.  It  should  be  noted,  however,  that  the  calculated  relative 
slopes  of  the  acoustic  and  erosive  burning  curves  do  not  correspond  to  the 
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experimentally  obtained  results. 

The  ability  of  the  acoustically  induced  turbulence  to  interact  with  the 
combustion  process  is  directly  related  to  the  height  of  the  acoustic  boundary-layer, 
Saai  which  is  inversely  related  to  the  frequency  of  the  acoustic  oscillations. 
Medvedev  and  Revyagin  reported  that  the  magnitude  of  the  acoustic  erosion  was  a 
weak  function  of  acoustic  frequency  and  suggested  it  may  be  due  to  the  decrease  in 
the  height  of  the  Stokes  layer.  Figure  31  shows  the  results  of  calculations 
performed  to  investigate  the  frequency  dependence  of  the  acoustic  erosion  for  the 
present  propellant  combustion  parameters.  The  calculations  were  performed  at 
f=1000  Hz  and  2000  Hz  with  an  equivalent  sand  roughness  height  of  ks=100pm. 
The  frequency  dependence  is  predicted  to  be  quite  strong  for  the  propellant 
thermophysical  parameters  of  this  investigation.  For  example,  at  Ue=60  m/s,  the 
acoustic  erosion  is  predicted  to  be  approximately  9  percent  at  f=1000  Hz  and  21 
percent  at  f=2000  Hz. 
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SECTION  4:  CONCLUSIONS 


An  analysis  of  transitional  and  turbulent  acoustic  boundary  layers  in  the 
presence  of  strong  injection  has  been  presented.  The  problem  was  approached  by 
analyzing  the  behavior  of  a  second-order  turbulence  model  rather  than  the 
traditional  Orr-Summerfeld  linearization.  An  approximate,  order-of-magnitude 
analysis  provides  a  simple  and  functionally  correct  estimate  of  transition  for  non- 
injected  acoustic  boundary-layers.  This  technique  may  prove  useful  for 
estimating  the  stability  characteristics  of  other  types  of  flows. 

The  computed  results  for  Stokes  layers  with  injection  indicate  a  substantial 
increase  in  the  acoustic  boundary  layer  thickness  for  strongly  injected  laminar  or 
turbulent  flows,  but  no  "blow-off'  condition  is  observed.  For  injected  fully 
turbulent  flows,  a  pronounced  Richardson  effect  is  obtained  and  is  accompanied 
by  very  large  maximum  levels  of  turbulence  (>  20%).  Both  axial  velocity  and 
turbulence  profiles  exhibit  a  convected  wave  train  shape.  The  turbulence 
development  was  also  found  to  lag  the  rectified  acoustic  velocity  by  a  few  degrees 
at  lower  frequencies. 

Calculated  transition  results  for  conventional  Stokes  layers  are  in 
agreement  with  the  approximate  analysis  and  data  trends,  although  at  the  higher 
acoustic  amplitudes  shown  by  the  scaling  relation,  additionsil  nonlinear  effects 
would  be  important.  The  effect  of  injection  is  to  (nonmonotonically)  decrease  the 
critical  acoustic  Mach  number  for  transition  by  up  to  a  factor  of  about  three  for 
smooth  surfaces.  This  effect  is  frequency  dependent,  but  is  expressible  in  terms  of 
the  acoustic  and  injection  Reynolds  numbers  for  the  problem.  The  results 
emphasize  the  importance  of  the  acoustic  and  injection  Reynolds  numbers  as 
principle  similtude  parameters  for  near-surface  phenomena. 

A  further  appreciable  reduction  in  critical  Mach  number  (or  Reynolds 
number)  is  predicted  for  injection  velocities  with  finite,  continuous  disturbance 
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levels.  This  effect  is  strong  enough  to  indicate  that  at  the  minimum  critical 
acoustic  Reynolds  number,  rather  modest  levels  of  acoustic  pressure  ratio  (~  1%) 
can  induce  significant  turbulence  levels  near  the  surface. 

The  potential  for  turbularization  of  a  near-surface  reaction  zone  in  ducted 
flows  such  as  rocket  chambers  therefore  theoretically  exists.  However,  whether 
this  mechanism  can  produce  instability  in  such  systems  depends  upon  several 
effects  not  considered  in  this  investigation.  For  example,  combustion  processes 
and  (acoustically)  nonlinear  behavior  such  as  thermoacoustic  streaming  could 
influence  the  nature  of  the  overall  response,  while  the  axial  mean  flow  could  bias 
the  transitional  characteristics  over  the  acoustic  velocity  and  pressure  cycles. 
Further  research  on  these  topics  is  in  progress. 

The  analysis  presented  in  this  work  has  attempted  to  identify  and  quantify 
some  aspects  of  the  role  of  acoustically  induced  turbulence  in  the  combustion 
response  of  solid  propellants.  Despite  several  simplifying  assumptions,  a 
complex  nonlinear  and  spatially  dependent  interaction  occurs  in  the  reaction 
zone  which  produces  an  analogously  complex  response  function.  Some  features 
of  the  response  have  long  been  hypothesized  (see,  for  example.  Ref  1.)  and  are 
presently  analytically  confirmed,  viz;  "full  wave"  rectification  relative  to  the 
harmonic  acoustic  velocity  outside  the  boundary  layer,  the  occurrence  of  a 
threshold  acoustic  amplitude  for  significant  response,  and  the  appearance  of  a 
time-mean  augmentation  ("D-C  shift")  of  the  propellant  regression  rate. 
Additionally,  the  magnitude  of  the  response  is  quite  sensitive  to  propellant 
"formulation  variables";  surface  roughness  and  gas  phase  heat  of  reaction 
(flame-zone  thickness)  being  those  assessed  in  this  study. 

The  results  also  indicate  differences  with  some  prevailing  conceptual 
aspects  of  velocity  coupled  response.  These  include: 


38 


1)  The  response  to  acoustic  velocity  oscillations  can  be  independent  of 
steady  state  erosive  burning,  as  demonstrated  by  the  neglect  of  an 
appreciable  axial  mean  flow  in  this  analysis.  Consequently, 
contemporary  treatments  of  velocity  response  functions  which  simply 
perturb  an  erosive  burning  response  would  omit  an  important 
additional  contribution  and  physical  mechanism. 

2)  For  nominal  combustion  parameters,  the  susceptibility  to  acoustic 
transition  and  propellant  response  increases  with  increasing  mean 
pressure. 

3)  The  appreciable  phase  shift  of  the  regression  rate  response  relative  to 
the  harmonic  external  axial  velocity  can  yield  a  significant 
'apparent'  pressure  response,  if  an  attempt  was  made  to  analyze  this 
with  conventional  linear  theory. 

4)  The  behavior  of  the  velocity  response,  even  at  moderate  acoustic 
pressure  amplitudes  of  a  few  percent,  is  highly  nonlinear  and 
possesses  significant  harmonic  content. 

The  present  results  have  been  obtained  utilizing  several  simplifying 
assumptions,  some  of  which  pose  significant  research  issues  in  their  own  right 
(for  example,  the  potential  for  direct  acoustic-turbulence  interaction,  or  the 
potential  for  extreme  combustion-turbulence  interaction  from  composite- 
propellent  diffusion  flamelets).  One  of  the  more  intriguing  questions  posed  by  the 
present  results  is  the  potential  influence  of  an  axial  mean  flow  interacting  with 
the  acoustic  boundary-layer  and  turbulence  fields.  In  this  case  a  substantial 
biasing  of  the  propellant  response  over  a  cycle  is  anticipated  (a  necessary 
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condition  for  instability).  Considering  the  sensitivity  and  nonlinearity  of  the 
acoustic  transition  process,  this  area  of  research  could  provide  special  insight 
into  the  coupling  mechanisms  of  producing  instability  in  actual  devices. 

The  results  presented  in  this  work  provide  an  attempt  at  evaluating  the 
combined  response,  due  to  simultaneous  acoustic  pressure  and  (turbulent) 
velocity  oscillations,  of  the  combustion  of  a  solid  propellant.  An  order  of 
magnitude  analysis  performed  on  the  mean  enthalpy  equation  demonstrated  that 
it  and  the  mean  species  equation  remain  similar  in  the  presence  of  acoustic 
oscillations  in  the  gas-phase  reaction  zone.  This  result  should  prove  useful  in 
future  investigations  of  similar  unsteady  problems  involving  premixed 
combustion,  since  chemical  species  transport  equations  do  not  have  to  be  solved 
computationally.  Despite  several  simplifying  assumptions,  the  acoustic  pressure 
and  the  turbulent  fluctuations  of  the  acoustic  boimdary-layer  interact  with  the 
reaction  zone  to  produce  an  intricate  response  of  the  propellant  combustion. 
Several  of  the  characteristics  of  acoustic  erosion  which  have  been  determined  or 
inferred  experimentally  using  T-bumers  are  presently  verified.  These  include: 

1)  The  maximum  of  acoustic  erosion  occurs  at  the  velocity  antinode 
(pressure  node)  and  the  minimum  occurs  at  the  ends  of  the  chamber 
(pressure  antinode).  The  calculations  did  not  show  the  negative 
acoustic  erosion  (pressure  response)  near  the  ends  of  the  chamber  as 
was  experimentally  reported.  The  authors  believe  this  discrepency  to 
be  the  result  of  a  simplified  modeling  of  the  propellant  combustion 
process. 

2)  Under  nominal  conditions,  the  threshold  velocity  of  acoustic  erosion 
is  less  than  half  the  threshold  of  steady-state  erosion. 
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3)  Acoustic  erosion  is  a  nearly  linear  function  of  the  amplitude  of  the 
acoustic  velocity  oscillation.  (However,  the  relative  slopes  of  the 
calculated  acoustic  and  steady-state  erosion  curves  did  not  agree  well 
with  the  experimentally  obtained  results,  again  a  possible 
consequence  of  the  combustion  model.) 

4)  The  magnitude  of  the  acoustic  erosion  is  a  function  of  the  frequency 
of  the  acoustic  oscillation. 

The  results  have  also  provided  the  variation  of  several  characteristics  of  the 
combined  response  as  a  function  of  axial  location  within  the  chamber.  The 
response  of  the  propellant  was  found  to  change  from  the  purely  sinusoidal 
response  at  the  chamber  ends  to  the  completely  rectified  response  containing  only 
even  harmonics  at  the  center  of  the  chamber.  Between  those  extremes,  the 
frequency  content  is  composed  of  the  fundamental  and  both  even  and  odd 
harmonics,  which  could  lead  to  the  excitation  of  other  modes  of  acoustic 
oscillation  within  the  chamber  of  an  actual  operating  system. 
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TECHNICAL  APPENDIX  A:  TABLE  OF  ASSUMPTIONS 


Summa^_^:Assum2tions_Concemin£jhe^Ga8^and^Conden8ed^hase8^^^^^^^^^_^^_^^ 

Gas  Phase  Assumptions 

a)  External  body  forces  and  the  coefficient  of  bulk  viscosity  are 
negligible. 

b)  Effects  of  radiation  are  negligible. 

c)  Species  diffusion  due  to  thermal  and  pressure  gradient  effects  is 
negligible,  and  all  binary  diffusion  coefficients  are  equal. 

d)  The  Lewis  niimber  is  unity. 

e)  Specific  heats  of  the  chemical  species  are  equal  and  independent  of 
temperature. 

f)  The  mean  flow  Mach  number,  <M>,  is  negligible. 

g)  The  acoustic  boundary  layer  thickness,  5a,  is  much  smaller  than  the 
duct  radius  or  transverse  dimension. 

h)  Combustion  of  a  homogeneous  reactant  mixture  proceeds  through  a 
single  step,  irreversible,  chemical  reaction. 

i)  The  gas  phase  species-average  molecular  weight  is  equal  to  the 
molecular  weight  of  the  mixture  at  equilibrium. 

h)  The  turbulent  enthalpy  correlations  above  the  flame  zone  are 
negligible. 

Condensed  Phase  Assumptions 

a)  The  condensed  phase  is  homogeneous. 

b)  The  condensed  phase  is  stationary  in  the  reference  coordinate 
system. 

c)  Material  properties  are  constant. 

d)  Species  diffusion  is  negligible. 
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e)  Radiation  effects  at  the  interface  and  in  the  condensed  phase  are 
negligible. 

£)  The  surface  regression  rate  is  small  compared  to  the  gas  phase 
injection  velocity. 

g)  Subsurface  reactions  are  not  present. 

h)  Surface  regression  mechanism  is  governed  by  the  Arrhenius 
P3rrolysis  law.   
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TECHNICAL  APPENDIX  B:  TABLE  OF  PROPELLANT  PARAMETERS 


Propellant  Thermochemical  Parameters 


Tf  =  2976K 
cp  =  1.92  X 103  J/kg-K 
Wav  =  25.8  kg/kg-mole 
TAg  =  2.0xl04K 

Bg  =  6.5  X  109  kgl-^  8-1  K-Pe 
<D  =  1.10 
Pg  =  0 

p  =  4.42  X  10-7  (T)0-65  kg/m/s/KO-65 
Pr  =  l 
T„  =  300  K 

Lg  =  —4.184  X  10®  J/kg  (-100  cal/gm) 

Prt  =  1700  kg/in3 

Cji  =  1.46  X 103 

Art  =  2.5  X  107  kg/m^-s 

TAn  =  1.0xl04K 

r  =  9.2  X  10-3  m/s 

T*=700K 

p*  =  68  atm 

n  =  0.55 
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Figure  1 .  Mean  flow  and  acoustic  profiles  at  pressure  node. 
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Figure  2.  Transition  to  turbulence  for  simple  acoustic  motion. 
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Figure  3. 


Figure  4, 
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Laminar  velocity  profiles  in  the  acoustic  boundary  layer 
without  injection. 
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Turbulent  velocity  profiles  in  the  acoustic  boundary  layer 
without  injection. 
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Figure  5. 


Figure  6. 
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Turbulence  intensity  profiles  in  the  acoustic  boundary  layer 
without  injection. 
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Turbulent  velocity  profiles  in  the  acoustic  boundary  layer  with 
injection. 
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Figure  7. 


Figure  8. 
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Turbulence  intensity  profiles  in  the  acoustic  boundary  layer 
with  injection. 
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Effect  of  injection  on  the  stability  of  the  acoustic  boundary 
layer. 
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Figure  9. 


Figure  10  (a-c). 


Effects  of  injection  and  injection  disturbance  level  on  the 
stability  of  the  acoustic  boundary  layer. 
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Time  variation  of  turbulence  development  for  the  injected 
acoustic  boundary  layer  (Oy  =  0,  y  =  5a). 
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Figvire  11  (a-c).  Time  variation  of  turbulence  development  for  the  injected 

acoustic  boundary  layer  (Oy  =  0.035,  y  =  8a). 
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Figure  12.  Steady  state  temperature  profiles. 
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Figure  13. 


Figure  14, 
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Comparison  of  temperature  and  turbulent  velocity  profiles  at 
various  times  within  a  cycle. 
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Comparison  of  corresponding  temperature  and  turbulence 
intensity  profiles  at  various  times  within  a  cycle. 
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Figure  15.  Critical  acoustic  Reynolds  number  required  for  transition. 
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Figure  16.  Variation  of  critical  pressure  ratio  required  for  transition  as  a 

function  of  mean  chamber  pressure. 
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Figure  17. 


Figure  18. 


Variation  of  propellant  response  as  a  function  of  time  for 
various  surface  roughness  heights. 
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Variation  of  propellant  response  as  a  funrtion  of  time  for 
endo-  and  exothermic  surface  decomposition  energies. 
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Figure  19. 


Figure  20. 
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Harmonic  content  of  the  propellant  responses  shovm  in  Fig.  16 
for  a  smooth  and  rough  wall  (kg  =  0  and  kg  =  150  nm , 
respectively) 
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Relative  phase  angle  of  the  responses  shown  in  Fig.  16. 
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Figure  21 . 


Figure  22. 
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Normalized  steady-state  temperature  profiles. 
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Comparison  of  normalized  temperature  perturbations  with 
and  without  the  similarity  assumption. 
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Figure  23. 


Figure  24. 
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Normalized  laminar  and  turbulent  gas-phase  temperature 
perturbations. 
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Normalized  laminar  and  turbulent  gas-phase  normal  velocity 
perturbations. 
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Figure  25. 
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Normalized  laminar  and  turbulent  condensed-phase 
temperature  perturbations. 
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Figure  26. 


Combined  response  of  the  condensed-phase  regression  rate  as 
a  function  of  axial  location. 
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Figure  27.  Harmonic  content  of  the  combined  response  versus  axial 

location  for  Hma  =  0.05,  f  =1000  Hz,  and  kg  =  150)im). 
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Figure  28.  Acoustic  erosion  vs  axial  location  within  the  chamber 

assuming  various  equivalent  sand  roughness  heights. 
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Figvire  29. 


Figure  30. 


Experimental  comparison  of  acoustic  and  steady-state  erosion 
(as  presented  in  Ref.  14.) 
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Comparison  of  calculated  acoustic  and  steady-state  erosion. 
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Figure  31 .  Comparison  of  calculated  acoustic  erosion  at  acoustic  oscillations  of 

1000  Hz  and  2000  Hz, 
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